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A systematic field theory is presented for charged systems. The one-loop level corresponds to the 
classical Debye-Hiickel (DH) theory, and exhibits the full hierarchy of multi-body correlations deter- 
mined by pair-distribution functions given by the screened DH potential. Higher-loop corrections 
can lead to attractive pair interactions between colloids in asymmetric ionic environments. The free 
energy follows as a loop- wise expansion in half-integer powers of the density; the resulting two-phase 
demixing region shows pronounced deviations from DH theory for strongly charged colloids. 



Since the early work of Debye and Hiickel (DH) it is 
known that electrostatic interactions in a mixture of pos- 
itively and negatively charged particles produce a net at- 
traction . This is due to charge screening: Each charge 
is (on average) predominantly surrounded by oppositely 
charged particles, which thus leads to an overall attrac- 
tion between the particles. The resulting DH free energy 
contribution has been theoretically demonstrated to lead 
to phase separation in the context of ionic fluids and 
colloidal mixtures || , in agreement with experimental Q 
and numerical J5| work on ionic fluids and experiments on 
colloidal mixtures Q . The exact nature and origin of the 
DH term has remained somewhat unclear, and several 
improvements have been devised based on series expan- 
sions and liquid-state theory [Q] , explicit incorporation of 
dipole pairs || , and density- functional theory || . 

A second intensely debated question concerns the pos- 
sible existence of attractive interactions between simi- 
larly charged objects in electrolyte solution ||. Experi- 
mentally, such an attraction has been seen for DNA ||] 
and strongly charged microspheres which are confined 
between charged plates [0. Clearly, the phase separa- 
tion observed for colloidal mixtures Q| is a priori not an 
indication for such an attractive interaction, because the 
dense phase is induced by attractions between oppositely 
charged particles, as becomes explicit within DH theory 
(also, see the discussion in f|). 

In this article we present a systematic field theory for 
charged systems, and calculate both the free energy and 
the effective interactions between charged particles im- 
mersed in an electrolyte solution. At the one-loop level, 
we recover the classical DH theory, the nature of which 
transpires in an especially lucid fashion within our frame- 
work: we find the full hierarchy of multi-body correla- 
tions to be present, with all pair-distribution functions 
given by the screened DH interaction. This means in 
specific that triplet correlations are already included at 
the DH level (in contrast to implicit assumptions in re- 
cent theories El ) , and that effective interactions between 
similarly charged particles are repulsive. At higher or- 
der in our theory (which corresponds to including multi- 



loop diagrams), non-trivial multi-body interactions ap- 
pear, and, consequently, the multibody correlations ac- 
quire contributions which cannot be described as super- 
position of pair correlations. Also, the effective pair in- 
teraction receives corrections which can be attractive if i) 
the electrolyte is asymmetric and consists of multivalent 
counterions and monovalent coions, or if ii) the colloidal 
charge is overcompensated by salt ions. The latter situa- 
tion is realized in experiments on charged microspheres, 
where a strong attraction is only found in the vicinity of 
a charged wall JlOj ] . The free energy of an ionic solution, 
expanded in the number of loops, follows to be a series in 
half-integer powers of the density, and thus constitutes 
a systematic low-density expansion The effects of 

higher-loop contributions on the demixing transition be- 
come increasingly important for highly charged colloids 
and lead to pronounced deviations from DH theory. 

To proceed, we consider the partition function of 
N charged, fixed test particles, immersed in a multi- 
component electrolyte solution with (in general) M dif- 
ferent types of ions, 
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where v(r) = £s/r is the Coulomb operator and the 
charge density operator p c is defined by 



M rij 



p c (r) ee W - R *) + E E ^ - r fc } ) ( 2 ) 

i=l j=l k=l 

with Qi and qj being the charges (in units of the el- 
ementary charge e) of the test particles and the ions, 
respectively. The length A is an arbitrary constant, 
and the Bjerrum length £b = e 2 /AireksT defines the 
length at which two unit charges interact with ther- 
mal energy fcgT. Electroneutrality of course requires 
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Noting that the inverse Coulomb operator can be explicitly written as v : (r) = — V 2 <5(r)/47r^s, after a Hubbard- 
Stratonovich transformation, the partition function is given by 



Z[{Rn}}= ||^exp| -^-|dr(V0) 2 -z^Q^(RO + X]^ lo g / 



N M r /■ , 
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where Z is the partition function of the inverse Coulomb operator, Z ~ detv, and the entropy of ideal mixing is 
S = — rtj ln(A Cj) with Cj = rij/V denoting the concentration of ion species j. Performing a cumulant expansion 
of (0) in powers of <f), we can rewrite the partition function as 



Z[{Rn}] = J |^ cxp | -~ J drdr>(r)^(r - r')<MO - i £ Q,^(R,) + + 5 j , 



(4) 



where i>dh is determined via the inverse operator equation (the so-called Dyson equation in field theory) 

vEk(r)=v-\T)+hS(T) (5) 

which is solved by the well-known DH interaction vdh(i") = £BC^ TK /r with the screening length k^ 1 defined by 
k 2 = 4tt£bI2- All anharmonic terms are contained in the non-local potential W, which is up to eighth order given by 

W[</>] = + ^ (f ~ S0 2 " 2 ) ~ l -^-(f -lOfi^j -^{^ ~10^ 2 + 3O0 2 " 3 ) 

We have introduced the generalized ionic strength I n , which is defined as /„ = ._ 1 g™Cj and can take both positive 

and negative values. In these equations (j> n denotes moments of the field, <fi n = f dr0 n (r) /V . The action is invariant 
with respect to a change of the gauge field <fr. We therefore set = 0. The linear term in <fi in Eq.(|J) can be removed 
by a shift of the fluctuating field 0, and the partition function then takes the form 



z \{Rn}] = exp 1 5 - ^Y^QiQjV-on^U - Rj)j /"^" exp {~^ /' 



drdr'0(i> D ^ (r - r')0(r') + V^[0] (7) 



where 0(r) = 0(r) — ij^i <9i u DH(r — Rj)- Up to this 
point, our calculations are (in principle) exact: keeping 
terms of all powers in W in ([?]) leads to a model equiva- 
lent to the original partition function ([j]). They are also 
systematic, in that keeping terms of higher and higher 
order of W should make the resulting theory a more 
and more faithful representation of the underlying phys- 
ical model. In fact, we will demonstrate that the lowest 
approximation, where the potential W and thus anhar- 
monic terms in are neglected altogether, is equivalent 
to the classical DH theory. In this case, it follows from 
that the dimensionless pair interaction C/ 2 between 
two test particles is just the DH potential, 

£/ 2 (Ri - R 2 ) = Q<Q 3 -«dh(Ri ~ Ha) (8) 

and the two-point correlation function is 52 (Ri — R2) °c 
e -c/(Ri-R 2 ) w ith a proportionality constant such that it 
is normalized Neglecting W, there are no multibody 
interactions between test particles in (Q), and higher- 
order correlation functions are therefore given by prod- 



ucts of the pair correlation function, 33 (Ri, R 2 , R3) oc 
02 (Ri - R 2 )ff2(R2 - R3)ff2(Ri - R3), and so on. This 
is the superposition principle, known as a postulate from 
liquid state theory; it is exactly obeyed in DH theory. To 
connect to liquid state theory, we note that from Eq.(||) 
one obtains by inversion the integral equation (which in 
fact holds for any pair interaction v) 

v(r) = t'DH (r) + h J dr'v DH (r')t'(r - r'), 

the field-theoretic version of the Ornstein-Zernicke equa- 
tion. The DH interaction is the exact solution of this 
integral equation Jl3| . The DH theory therefore contains 
correlations of all orders, and there is no need to explicitly 
add higher-order correlations (compare ||). Improve- 
ments can only come from adding non-trivial higher-body 
effective interactions, i.e., from violations of the superpo- 
sition principle. This is the effect of the potential W, as 
we will demonstrate in the following. 
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Expanding W[<p] in the exponential of Eq. ([?]), the first correction comes from the cubic term, 

Z[{R N }] oc exp \ S- ^^2QiQjV DH (Ri -R-i) - QiQjQfcQ3(Rj,Rj)Rfc) 

with the three-point vertex given by f2 3 (Ri, R2, R3) = / dr VDn(r — Ri)udh(i" — R2)vdh( f — R-3)- 



(9) 



The summation over (i,j,k) in (|]) is unrestricted. By 
considering the case where two of the three indices are 
equal, one obtains a correction to the pair interaction, 
which reads 



The DH partition function is Z% ~ detuDH, whereas 
higher-order correlations are contained in W. The expec- 
tation value in ([l2]) is evaluated with the DH propagator. 
We first evaluate the DH free energy, 



AC/ 2 (R) 
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where the function H (which is positive) is determined 



by the integral S(x) 
27r{e- a: (ln3-r[0,a;]) 
behavior of H is 
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where the momentum integral goes over a cube of length 
2tt J a. Since the integrand is isotropic we distort the in- 
tegration volume to a sphere and obtain after a straight- 
forward integration 
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and thus shows the same asymptotic behavior as the DH 
repulsion for large separations. When do we expect at- 
tractive interactions between likely charged particles? In 
other words, when does the prefactor in Eq.([l0|) become 
negative? If we assume the test particles to be posi- 
tively charged, the condition for attraction is that the 
third-order ionic strength I3 is negative. Assuming a 
homogeneous solution of positive macroions with charge 
Z, concentration c, and counter ions of valency z, the 



third-order ionic strength is I3 = cZ{Z 



and clearly 



always positive: similarly charged particles at finite con- 
centration do not attract each other, in agreement with 
experiments. On the other hand, considering two single 
charged macroions, it is easy to see that ^3 is negative if 
the salt solution is asymmetric: for positive macroions, 
attraction is therefore possible i) if one has an electrolyte 
consisting of negative ions with a higher valency than 
the positive ions or ii) if there are more negative than 
positive ions in the local environment (as in experiments 
between two charged walls Eol). Comparing the strength 
of the DH repulsion Eq.(||) and the attraction Eq.(10) at 
large separation we find the attraction to dominate for 
c£%Z 2 > 9(TO 2 + l)/(m 3 -l) 2 7rln 2 3forthe case of a 7Ti : 1 
electrolyte jiS. Experimentally, it is well-known that 
asymmetric salts like Calciumchloride induce the precip- 
itation of negatively charged macroions and negatively 
charged polymers |J. 

The phase behavior of charged colloidal mixtures fol- 
lows from the free energy with all particle coordinates 
integrated over, 



T = - hxZ= -S -In 
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In the limit a — > one obtains the well-known result 
/dh — — a 3 K 3 /127r (plus corrections which scale linearly 



in k and thus correspond to an unimportant shift in 
the chemical potential). Our corrections as a function 
of the cut-off take in an approximate fashion the fi- 
nite ion-sizes into account. The free energy contribution 
A/ = -^ln(e w M) is, using Eq.(|), given by 
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The expectation value (</> 2 ) is 
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and the generalized susceptibility \ n is defined as 



Xn 



dr((/> <M n = 47r^(77K)"- 3 r[3 - 77, ann] 



Naive scaling predicts the generalized susceptibilities in 



A/ to scale like Xr. 



2 (n 3)/2 ag a mnc tj on f t ne j on 



density c. Since the dominant terms in A/ scale as I^Xn, 
one would thus obtain a systematic free-energy expansion 
in half-integer powers of the density, starting with the 
DH term /dh, which asymptotically scales as c 3 / 2 jn|]. In 
practice, the integrals in \n diverge in the ultraviolet and 
thus depend on the ion radius a in a crucial way, which 
leads to changes from the naive scaling picture. We now 
present results for the simplified case of colloids of charge 
Z and concentration c with counterions of valency z. In- 
troducing the energy scale e = zZ£s/a and the total ion 
volume fraction c = a 3 c(l + Z/z) the DH theory (which 
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amounts to neglecting the term A/) predicts a critical 
point at e ~ 5.63 and c ~ 0.0418, independent of the col- 
loid charge Z . Including higher-order terms contained in 
A/ the critical interaction strength e and volume fraction 
c depend on Z/z, as shown in Fig. 1. The limiting values 
for Z/z = 1 are e ~ 5.61 and c ~ 0.0416 and are thus 
very similar to the DH case [|l5| . For Z/z > 1 the devia- 
tions from DH theory are pronounced. In Fig. 2 we show 
coexistence curves for Z/z — 1, 2, and 10 (solid, broken, 
and dotted lines, respectively). The coexistence curve as 
predicted by DH theory would be indistinguishable from 
the solid line. 

In summary, we introduced a systematic field theory to 
describe charged colloidal suspensions, including fluctu- 
ations of and correlations between charged particles. At- 
tractive colloidal interactions are predicted for asymmet- 
ric electrolyte solutions, as realized for multivalent m : 1 
salt solutions and in the neighborhood of charged walls. 
A critical point of demixing occurs even in the absence 
of attractive colloidal interactions, whose location in the 
temperature-density plane depends strongly on the col- 
loidal charge. The hard-core repulsion between colloids 
and ions has been incorporated by a small-distance cut- 
off, similar to recent lattice theories [El| . This is a rather 
poor description of colloidal systems, because here a large 
size difference between colloids and ions exists. We hope 
to treat size asymmetries more accurately in the future 
using perturbative treatments of the hard-core interac- 
tions between all particle pairs. 
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FIG. 1. Critical interaction strength e = zZls/a (solid 
line) and critical volume fraction c = a 3 c(l + Z/z) (broken 
line) as a function of the colloid/counterion charge ratio Z/z. 

FIG. 2. Two- phase coexistence regions for Z/z — 1 (solid 
line), 2 (broken line), and 10 (dotted line) as a function of 
volume fraction c and rescaled temperature z/Ze. 
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Fig. 2 




